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Abstract — The paraxial lens data were determined for a series of 
symmetrical magnetic lenses of equal lens diameter but variable air 
gap width for a wide range of lens excitati ons using the three- 
dimensional electrodynamic computer code MAFIA The results are 
compared with a similar study done by Liebman and Grad wherein the 
fieto distributions within the lenses were measured experimentally with 
a resistance network analogue Using these fields the lens dataware 
obtained through numerical trajectory tracing. The utility of using 
MAFIA, instead of experimental methods for lens design is shown by 
the excellent agreement of the simulated results compared to 
experiment Also demonstrated is the capability of using MAFIA to 
investigate aberration sources such as higher order off-axis magnetic 
field and space-charge effects. 

L Introduction 

Paraxial lens characteristics were simulated for a series 
of common symmetrical magnetic lenses using the computer code 
MAFIA (MAxwelTs Equations solved by the Finite Integration 
Algorithm). The results are compared with the experimentally 
obtained results published in [1] with excellent agreement 
Thus, the accuracy of MAFIA indicates that time-consuming and 
costly experimental test measurements in lens design can be 
greatly reduced. 

MAFIA is a powerful, three-dimensional electromag- 
netic code written in FORTRAN 77. It is used for the computer- 
aided design of fully 3D and 2D electromagnetic devices, magnets, 
RF cavities, waveguides, antennas, etc. The Finite Integration 
Technique (FIT) algorithm produces a set of finite-difference 
matrix equations for the electric and magnetic field vectors in the 
structure under study. The solution of these equations yields 
static, fiequency-domain or time-domain solutions of Maxwell's 
equations [2], [3]. The code includes nine interrelated modules. 
The lens data of this repent were calculated using the M (mesh 
generator), S (static solver), TS2 (two-dimensional (2D) 
particle-in-cell (PIC) solver), and P (postprocessor) modules of 
MAFIA. The remaining modules include an eigenmode solver, 
3D PIC solver, 2D and 3D time domain solvers, and a frequency- 
domain eddy-current solver. 

To test the utility of MAFIA, a study was conducted 
simil ar to [1] wherein a series of magnetic lenses was investi- 
gated. Since [1] indicates that the main geometrical parameter of 
a symmetrical lens is the ratio S/D (see Fig. 1), the dependence 
of the paraxial ima ging properties of a series of lenses was 


computationally investigated where the coil diameter D was held 
constant and the air gap width S varied. The field distributions 
within the lenses calculated with the magnetostatic solver of 
MAFIA are compared to fields obtained experimentally with a 
resistance network analogue. The trajectories of electrons emitted 
from a cathode surface, accelerated through a constant potential 
and focused by the simulated magnetic electron lenses were 
calculated with the 2D PIC module of MAFIA and compared to 
paths calculated in [1] by numerical trajectory tracing. Aberration 
sources such as off-axis fields and space-charge effects are also 
computationally examined, demonstrating how these very 
important parameters in lens design are simulated using the 
MAFIA code. 

IL Analysis 

The 2D MAFIA plot (generated using the M module) for 
the magnetic electron lens (axially symmetric about the zz* axis) 
is shown for the r-z plane in Fig. 1 . A cathode and anode are 
placed at opposite axial ends ofa stainless steel tunnel. The static 
solver module S has the ability to solve electrostatic and 
magnetostatic problems. The electrostatic solver was used to 
establish a constant potential V f between the cathode and anode 
and calculate the resulting electric field. The magnetostatic 
solver of MAFIA was used to calculate the magnetic field 
associated with the lens. A soft magnetic material was used for 
the magnetic circuit and a current density was defined across the 
cross-sectional winding area representing the current windings 
with ampere-tum product NI (See Fig. 1 ). 

The MAFIA module TS2 computes the time integration 
of electromagnetic fields selfconsistently with the time integration 
of the equations of modem of charged particles that move under 
the influence of those fields. Since the fields caused by those 
moving charges are also taken into account, effects like space 
charge and magnetic forces between particles are fully simulated 

The electric fields calculated through the electrostatic 
solver and the magnetic fields calculated through the mag- 
netostatic solver were loaded into the TS2 simulation l A particle 
with the physical properties of an electron was emitted from the 
cathode surface. Such properties as particle position and dynamic 
electric field were monitored These monitored values can be 
displayed or subsequent calculations can be done on the fields 
using the postprocessor module P. 




Fig. 1 2D MAFIA plot in cylindrical coordinates for common 
electron lens s ymm e tri cal about the zz* axis 

m. Results 

A Static Calculations 

Thedimenskmlessexe^ P is presented in 

Fig. 2 of [1] as a function of the relative pole piece gap width 
S/D, where 
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where V, is the relativistically corrected accelerating voltage [4], 
R is the coil radius in m, and B c is the maximum axial field 
strength in Tesla. Fig. 2 compares the experimental B. calculated 
from (2) with the field strength calculated with MAFIA as a 
function of S/D. The agreement is very good, die small 
discrepancy probably due to the lack of information about the 
(B.H) curve of the iron used for the experimental lens pole pieces, 
pole piece shape and exact lens dimensions. Fig. 3 compares die 
experimental halfwidth to the values calculated from the MAFIA 
field distribution. The halfwidth a is defined as the z value for 
which the field strength has Men to BJ2. Both Figs. 2 and 3 
show that the on-axis field distribution for the series of coils is 
indeed accurately simulated with MAFIA. 

B. Particle-ln-Cell Calculations 

The particle positions were monitored through time for 
an electron emitted into the focusing configuration of Fig. 1 at a 
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Gap width S/D 


Fig. 2 MAFIA and Experimental peak magnetic flux density on the axis as 
a function of fee gap width S. 



Gap width S/D 


Fig.3 MAFIA and Experimental half-widths a as a function of the gap width 

S. 

radius r = 0. 1R, parallel to the zz'-axis. Fig. 4 shows the focal 
points (defined as die intersection of die trajectory with the 
"optical" axis), calculated by numerical trajectory tracing in [1] 
versus those simulated with MAFIA as a function of the excitation 
parameter k* of (2). The fiscal points of MAFIA are very close 
to those cal c ula te d with numerical methods. It is seen that as the 
accelerating potential is increased (k* is decreased), the particle's 
focal point moves closer to the electron source (z/R—2.5). The 
largest discrepancy occurs for excitation k*= 100. In this study, 
this corresponds to an accelerating potential of about 100 V. At 
such a low value of voltage, die magnetic effects of the lens coil on 
the trajectories is more prominent than the electrical effects. The 
discrepancy that was seen in Figs. 2 and 3 between the experimen- 
tal magnetic flux density and die MAFIA calculated values will 
have the largest effect for this case. Furthermore, inaccuracies in 
the numerical trajectory tracing will also be most evident 







Fig.5 shows the effects cf spherical aberrations occurring 
when particles are emitted at different radii. From the MAFIA 
trajectory plots, it is seen that the particles emitted at a larger 
value of r/R are focused closer to the source. This effect can be 
explained by the high-order terms of the off-axis magnetic fields 
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MAFIA focal points and those obtained through numerical trajectory 
tracing usi ng experimentally obtained field distributions for m a gn etic 
kns with gap width S/D= 1.0 as a function of the excitation parameter 

k*. 


forces present in a high -density electron source, a particle bunch 
consisting of 2000 particles with current density 100 A/cm 2 was 
emitted through the same magnetic lens configuration and its path 
monitored. Fig. 6 shows the MAFIA plot of the particle bunch and 
foe anowplot of the dynamic electric field, the size of the arrows 
proportional to foe magnitude of foe field. At this instant in time 
it can be seen that the particles have diverged from their initial 
radial positions and acquired various radial and axial velocities. 




Axial position z/R 


Fig. 5 MAFIA trajectory pads of electrons emitted separately at various radii 

through imagjttfic lens vrith gap width S/I>=1.0 and excitation 1^5.0 
illustrating spherical aberration effects. 

which are often excluded in basic electron-optical theory [4], [5], 
but are included in the MAFIA simulations. 

hi a large number of electron lens focusing applications, 
a high density electron source is used With such a source, 
collisions between electrons, together with the Coulomb forces of 
repulsion tend to diverge the electrons. Ifthereisno neutraliza- 
tion of this electron space charge by ions, the radial spread will 
introduce aberrations [4]. Thus, in this case the trajectory paths 
will not follow basic electron-optical theory as closely. To 
demonstrate the ability of MAFIA to simulate the space charge 


Fig. 6 MAFIA plot of high-density electron bunch and the associated 

dynamic electric field through a magnetic lens with gap width S/D“=1.0 
and excitation k*~5.0. 


IV. Conclusions 

A series of common symmetrical magnetic electron 
lenses was simulated with the computer code, MAFIA. The 
m magnetic flux density was calculated with the magnetostatic solver 
of MAFIA for the series keeping the lens diameter constant and 
varying the air gap width for a large range of excitations. The 
MAFIA results are compared to experimental values calculated 
with a resistance network analogue with good agreement The 
qrmilfltfri magnetic field was subsequently loaded into foe 2D PIC 
module of MAFIA and foe trajectory paths through a constant 
accelerating field calculated. When compared with the focal 
lengths obtained from numerical trajectory tracing using the 
experimentally obtained field values, the focal lengths obtained 
with MAFIA showed good agreement with the largest discrepancy 
occurring at a potential of about 100 V. This deviation could be 
due to foe difference in the simulated magnetic flux density of the 
coil compared with experiment or errors in the numerical 
trajectory tracing of [1] itself due to simplifications involved in 
the calculations. With more information about the specific lens 
parameters such as dimensions and magnetic material, it is 
believed that the accuracy could be improved. 
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